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Abstract 

In heavy-ion collisions at high energies, the quantum chromodynamics (QCD) pre- 
dicts the production of the deconfined quark-gluon plasma (QGP) state. Quarkonia 
(cc or bb bound states) are a useful means to probe QGP and to investigate the 
behavior of QCD under the high parton-density environment. Up to now, the large 
hadron collider (LHC) at CERN provided two runs for PbPb collisions at a/snn = 
2.76 TeV in the years 2010 and 2011. The ALICE, ATLAS, and CMS experiments at 
LHC have analyzed the yields and spectra of the // ip and Y families. In this article, 
we review particularly the recent charmonium results in PbPb collisions at LHC from 
the 2010 run. 



Presented at Charm2012: The 5th International Workshop on Charm Physics, University 
of Hawai'i, Honolulu, U.S.A., May 14-17, 2012 



1 



1 Introduction 

QCD, the theory of strong interaction, predicts the deconfined quark-gluon plasma (QGP) state 
of partonic matter in high energy heavy-ion collisions. In particular, the heavy-ion collisions at 
colliders, such as LHC, create QGP in high temperature and small net baryon density, which 
resembles the primordial matter of the universe that existed about a few jis after the Big Bang. 
The temperature of such QGP should be at least 180 MeV, the transition temperature, with the 
energy density larger than 1 GeV/fm 3 H). 

Among several experimental signatures for QGP, the suppression of quarkonium states is one 
of the most prominent observables Q. In the original picture, the number of quarkonia, created 
in primary nucleon-nucleon (NN) collisions, should be significantly reduced by the dissocia- 
tion of QQ (Q = c or b) bound states during the deconfined phase of quarks and gluons due 
to the color screening of the long-range QCD confining potential. Various quarkonium states 
are expected to be sequentially dissolved due to the color screening effect, depending on their 
binding energies and the temperature of the surrounding matter (HUH. As a result, quarko- 
nia can play the role of a thermometer for the temperature reached by heavy-ion collisions 
and provide valuable information on the detailed properties of QCD in the high-density phase 
0. Table, [I] summarizes the properties of ] /\p(lS), XcO-P) an d ip f (2S) obtained by the non- 
relativistic potential model [3, 4J. Since the binding energies and radii are different for various 
cc states, a sequential suppression pattern similar to Fig. [I] is expected [3J. 



State 


j Wis) 


Xc(lP) 


t//(2S) 


Mass (GeV/c 2 ) 


3.10 


3.53 


3.68 


Binding energy (GeV) 


0.64 


0.20 


0.05 


Radius (fm) 


0.25 


0.36 


0.45 



Table 1: Masses, binding energies, and radii of various charmonium states from non-relativistic 
potential model OH). 
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Figure 1: Sequential suppression pattern of the cc states expected in high-energy heavy-ion 
collisions |3|. 

However, there are some complications in interpreting the quarkonium production data in not 
only heavy-ion collisions, but also elementary pp collisions. First of all, at present, there is 
no satisfactory model that can successfully describe both the cross section and the polariza- 
tion of quarkonia even in pp collisions. Secondly, no detailed information on the cold nuclear 
matter effects, like shadowing and gluon saturation in initial states, are available with enough 
precision [6J. Thirdly, a similar J/xp suppression at RHIC (-^s^iv = 200 GeV) and SPS (y/s^j 
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Figure 2: Comparison of the nuclear modification factor of J/xp as a function of the event cen- 
trality between the SPS (t/s^ = 17 GeV) and RHIC (t/s^ = 200 GeV) data I7HTOI. At RHIC, 
the data at midrapidity (open circles) and forward rapidity (closed circles) are displayed. The 
nuclear modification factor Raa is the particle yield in nuclear (AA) collisions relative to that 
in pp collisions after proper scaling by the number of binary collisions. (The definition of Raa 
will be described in detail in Sec.[3j) N part is the number of participating nucleons determined 
by the Glauber model so that a larger N part means more central collisions with a smaller impact 
parameter. 

= 17 GeV) and a larger suppression at forward rapidity than at midrapidity in AuAu colli- 
sions at RHIC still remain as puzzles (Fig. [2]) H7l-H0ll. Finally, the model calculations are not 
converging to the measured J/xp productions in nuclear collisions. For example, the statistical 
coalescence (or recombination or regeneration) model of uncorrected c and c quarks predicts 
a slight increase, but the holographic QCD (AdS/CFT) calculation with J/xps embedded in a 
hydrodynamic model predicts a large decrease of the J/xp yield with respect to the pp data 
for the transverse momentum (pj) larger than ~5 GeV/c in CuCu collisions at RHIC |TTl[T2| . 
As a result, the quarkonium production data at much higher LHC energies must be critical to 
resolve those questions. 

So far, LHC has provided PbPb collisions at t/snn = 2.76 TeV in two periods. During the first 
heavy-ion run in 2010, LHC delivered about 7 }ib~ l of integrated luminosity. In the follow- 
ing year 2011, thanks to the great improvement in the accelerator performance, LHC delivered 
about a factor of twenty more (about 150 }ib~ l ) PbPb events during a similar data taking period 
as 2010 (Fig. [3]). With much higher statistics of the new data sample accumulated in 2011, the 
detailed analyses for xp f (2S) and higher Y states become possible, but in this article, we present 
the experimental data only from the 2010 run. Note that LHC has also delivered pp collisions 
at the same energy as PbPb in March 2011, and the experiments used them as references when 
they analyzed the PbPb data. For example, CMS has accumulated about 225 }ib~ l of the in- 
tegrated luminosity of pp (C vv ), which is roughly the same number of binary collision scaled 
PbPb data recorded in 2010. 

This paper summarizes the charmonium production in PbPb collisions at y/sjsjN = 2.76 TeV 
from the ALICE |13|, ATLAS d, and CMS ED collaborations at LHC. The results on non- 
prompt J/xp from B-meson decay lH5ll and Y states lH6ll from CMS are also briefly presented. 
The detailed descriptions on the ALICE, ATLAS, and CMS detector systems can be found in 
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Figure 3: Integrated luminosities delivered by LHC (top) and recorded by CMS (middle) in 
2011 as a function of Day /Month. For comparison, the integrated luminosities for the first 
heavy-ion run in 2010 are also displayed in the bottom. About a factor of twenty enhancement 
was achieved in the integrated luminosities at LHC. 



2 J/ip at High-p r from ATLAS 

Soon after the first heavy-ion run at LHC, the ATLAS collaboration has deliberately published 
the centrality dependence of the J/xp yield [14J. Figure [4] displays the invariant mass distribu- 
tion of pairs in the // xp region for the most central class (0 - 10%) of PbPb collisions. Note 
that about 80% of the J/xps reconstructed by ATLAS are in pj > 6.5 GeV/c. 

Since the reference pp data were not available at the time of this analysis, ATLAS estimated the 
relative // xp yield in each centrality bin 'cent' by normalizing the yield in the most peripheral 
class (40 - 80%), N^ h (J / xp) / N^ 0% (J / xp) . Following this, the central-to-peripheral ratio, R cp , 




\i + jr invariant mass [GeV] 



Figure 4: Invariant mass distribution of ji^ji pairs in the // xp region for the most central bin 
(0 - 10%) in PbPb collisions at ^/s^ = 2.76 TeV from ATLAS EH. 
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3 Jlip at High-pr from CMS 
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Figure 5: R cp of J/xp as a function of centrality in PbPb collisions at t/snn = 2.76 TeV from 
ATLAS |14|. The vertical bars represent the statistical errors while the grey boxes represent the 
combined systematic errors. 

was obtained by scaling the relative yield by the ratio of the mean number of binary collisions, 
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Figure [5] summarizes R Cjt7 as a function of centrality. It clearly demonstrates that the yield of 
J/xp is suppressed for more central collisions with a smaller impact parameter. Note that the 
uncertainty displayed with the most peripheral point, which is equal to unity by construction, 
is global to all points. 



3 J Ixp at High-p r from CMS 

The geometrical acceptance of the CMS detector is similar to that of the ATLAS detector for 
muons. The low pj limits of CMS for the accepted muon pairs are 6.5 GeV/ c for the midrapid- 
ity bin (| y | < 1.2) and 3 GeV/c for the most forward rapidity bin (1.6 < |y| < 2.4) |15|. CMS 
used the pp data as references and, instead of R cp , analyzed the nuclear modification factor 
Raa'. 

_ ^VV ^PbPb(//^) e vv £2) 

TaaNmb N pp {J/xp) epbpb 7 

where Taa is the nuclear overlap function estimated by the Glauber model, Nmb is the count 
of equivalent minimum bias events in PbPb, and e^/epbPb is the multiplicity dependent ratio 
of the efficiencies in pp and PbPb collisions for trigger and reconstruction. The ratio, e^/epbPb/ 
was determined by the Monte-Carlo simulation that embedded a PYTHIA signal event l20ll 
to a HYDJET background event |2T|. In the simulation, e^/epbPb was estimated to be about 
1.17 for the most central bin. By definition, deviation of Eq.|2]from unity indicates the medium 
effect on the J/xp production in heavy-ion collisions. 

In addition, CMS classified the measured J/xps into the prompt component associated with the 
primaty vertex and the non-prompt component associated with a secondary vertex of y+yr 
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Figure 6: Invariant mass (left) and pseudo-proper decay length (right) distributions of 
pairs for | y | < 2.4 in minimum bias PbPb collisions at y / S]\nv = 2.76 TeV from CMS IH5I . In the 
left panel, the solid line represents the sum of the 'crystal ball' function for the signal and an 
exponential function for the background. In both panels, the dotted and dashed lines represent 
the backgound and the non-prompt J /xp component, respectively. The solid lines are the total 
fit functions including the prompt J/xp component. 



|[T5ll . The prompt and non-prompt components have different origins: the former is from either 
direct production or feed down of the higher states, such as xp'(2S) and x c (lP), and the latter 
is from the B-meson decay. Because B mesons fly finite pathlength before they decay, the non- 
prompt J/xps can be separated from the prompt ones by the pseudo-proper decay length: 

/ - t m]l ^ n\ 
ij/tp - Lx y~^> 

where L xy is the distance between the primary and the secondary vertices in the transverse 
plane and tnj/y is the rest mass of J/tp. CMS has fitted the invariant mass and the pseudo- 
proper decay length distributions simultaneously to determine the yields of the prompt and 
non-prompt components. Figure [6] shows the invariant mass and the projected Ij/^ distribu- 
tions of fi + }t~ pairs in PbPb collisions at = 2.76 TeV [15J. Note that the resolution of the // ip 
peak is ~34 MeV/ c 2 , which is compatible to that in pp collisions. 

The fractions of B-meson decays in inclusive J/xps in pp as well as PbPb collisions at ^/sjvn = 
2.76 TeV are presented in Fig. [7] as a function of pj, integrated over centrality. For comparison, 
Fig. [7] also includes the data for b fraction measured by CDF in pp at y/s = 1.96 TeV l22l and 
CMS in pp at y/s = 7 TeV l23l . Within uncertainties, a good agreement between the present data 
and earlier results can be found. 

The nuclear modification factor Raa of prompt // xp is displayed in Fig. [8] as a function of the 
participating nucleons N part |15|. Raa in Fig- [8] is integrated over 6.5 < pr < 30.0 GeV/c in 
|y| < 2.4, and a large dependence on the collision centrality is observed. Raa of prompt // xp 
is about 0.6 for 50 - 100% centrality bin, and decreases to 0.2 ± 0.03(stat.) ± O.Ol(syst) ± 
0.01 (global), which indicates a factor of five suppression relative to pp, for the most central bin 
(0 - 10%). 

Figure [8] compares the current CMS results with the STAR preliminary data at RHIC {pj > 5 
GeV/c in |y| < 1.0) measured in AuAu collisions at y/smt = 200 GeV |24]|. This comparison 



6 



3 J/xp at High-p T from CMS 
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Figure 7: Fraction of the B-decay contribution to the inclusive J/xp yield as a function of pj 
in pp and PbPb collisions at ^/ijw - 2.76 TeV from CMS [15\. For comparison, the b fraction 
measured by CDF in pp collisions at y/s = 1.96 TeV 023 and by CMS in pp collisions at y/s = 
7 TeV II23II are also displayed. Points are plotted at the measured average values of pj. The 
statistical and systematic uncertainties are shown by bars and boxes, respectively 



reveals a larger suppression at higher beam energy at LHC for all centrality bins. Furthermore, 
CMS finds that the J/xp suppression apparently exists even in the most peripheral bin, which 
could not be observed in the R cp analysis done by ATLAS. 

Figure [9] shows the pr and rapidity dependences of Raa for prompt J/xp in PbPb collisions at 
y/s^ = 2.76 TeV from CMS, integrated over - 100% of centrality [HQ. Note that the data in 
Fig. [9] is integrated over |y| < 2.4 for the pj dependence in the left panel and over 6.5 < pr < 
30.0 GeV/ c for the y dependence in the right panel. Suppression by a factor of three is observed 
for two pr bins when integrated over centrality and rapidity without particular observable pj 
dependence. For the pj dependence of Raa the CMS data show a larger suppression at LHC 
energy than RHIC energy for pj > 5 GeV/ c Il25ll . 

In the rapidity dependence of Fig. [9} the CMS data show that the yield of prompt J/xp is less 
suppressed at forward rapidity than midrapidity in pj > 6.5 GeV/ c. This trend is opposite 
to what PHENIX has observed at the lower pj region (< 4 GeV/c) as shown in Fig. [2] [25 j. 
The difference between the CMS and the PHENIX data may come from (anti-)shadowing for 
different x ranges covered and/or different contributions by the regenerated // xps. However, 
further investigation is needed for the quantitative interpretation of the experimental data. 

Figure [10] displays the nuclear modification factor Raa of non-prompt J/xp as a function of 
Npart in PbPb collisions at y^jw = 2.76 TeV IH5I . It shows that Raa of non-prompt J/xp is 
0.37 ± 0.08(stat.) ± 0.02(syst.) ± 0.02(global) for the most central bin (0 - 20%) in PbPb colli- 
sions, which implies the suppression by a factor of ~2.7 relative to the pp data. This is the first 
data hinting the b-quark energy loss in medium. No centrality dependence is noticed with the 
present wide bin size, especially, in the peripheral bin of 20 - 100%. The data with finer central- 
ity bins are expected from the ongoing analysis using higher statistics data sample obtained in 
2011. 
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Figure 8: Nuclear modification factor Raa of prompt J/ip as a function of N par t in PbPb colli- 
sions at y/sjm = 2.76 TeV from CMS [15J. The measured // ips are for pj from 6.5 to 30 GeV/ c in 
|y| < 2.4. The grey box at R^ = 1 represents a global scale uncertainty of 6% originating from 
the integrated luminosity of the pp data sample. The statistical and systematic uncertainties 
are represented by bars and boxes, respectively. The CMS data are compared to the STAR data 
(stars) for p T > 5 GeV/c and |y | < 1.0 in AuAu collisions at y^NN = 200 GeV l24ll . 
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Figure 9: Nuclear modification factor Raa of prompt // xp, integrated over centrality, as func- 
tions of p T (left) and y (right) in PbPb collisions at v^nn = 2.76 TeV from CMS HH|. R A a is 
integrated over \y\ < 2.4 for the pj dependence in the left plot and over 6.5 < pr < 30.0 GeV/ c 
for the y dependence in the right plot under the minimum-bias centrality condition. The grey 
boxes at Raa = 1 represent global scale uncertainties of 6% originating from the integrated lu- 
minosity of the pp data sample. The statistical and systematic uncertainties are shown by bars 
and boxes, respectively. The PHENIX and STAR data obtained at RHIC are also included for 
comparison [24112511 . 
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4 J/tp at Low-p T from ALICE 
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Figure 10: Nuclear modification factor Raa of non-prompt J/xp as a function of N par t in PbPb 
collisions at y/sjm - 2.76 TeV from CMS |15|. The pj range of J/xp is from 6.5 to 30.0 GeV/c in 
|y| < 2.4. The grey box at Raa = 1 represents a global scale uncertainty of 6% originating from 
the integrated luminosity of the pp data sample. The statistical and systematic uncertainties are 
shown by bars and boxes, respectively. For comparison, the Raa data of prompt J/xp obtained 
by CMS are also included after rebinned as non-prompt J/xp in N par t- 



4 J/tp at Low-/7 r from ALICE 

Because the muon spectrometer is positioned in a very forward region, ALICE can measure 
single muons down to pj - 0.5 GeV/c and J/xps down to pj - GeV/c in 2.5 < y < 4 UTTll . 
Therefore, the geometrical acceptance of ALICE is quite complimentary to those of ATLAS and 
CMS. 



Figure 11 shows the invariant mass distributions of pairs from the same events and the 

mixed events in the J/xp region I26ll27| . Here, the invariant mass distributions from the same 
events include the J/xp signal as well as the backgrounds whereas those from the mixed events 



provide only the shape of the backgrounds. In particular, Fig. 11 compares the invariant mass 
distributions for the semi-central class (50 - 80%) with that for the most central class (0 - 10%) in 
PbPb collisions at ^/s^ = 2.76 TeV from ALICE. As expected, the signal-to-background ratio 
is significantly reduced for more central collisions. 

Figure |12| shows the nuclear modification factor Raa f° r inclusive J/xp in 2.5 < y < 4 and pj > 
GeV/c measured by ALICE in PbPb collisions at ^/sjviv = 2.76 TeV [13J. The centrality inte- 
grated Raa is 0.545 ± 0.032(stat.) ± 0.084(syst.), which clearly indicates the // xp suppression 
relative to the pp data. The contribution from B-meson decays to the inclusive // xp yield in the 
ALICE acceptance has been estimated using the LHCb measurements in pp collisions at y/s = 
7 TeV. The LHCb collaboration estimated the b fraction of the inclusive J/xp yield to be about 
10% in pp collisions at 7 TeV. Assuming the same fraction in PbPb at 2.76 TeV, ALICE has found 
that the difference in the Raa results between the prompt and inclusive J/xp does not exceed 
11%. 



For comparison, Fig. 12 also displays Raa of J/xp from PHENIX in AuAu collisions at ^Sjviv 
= 200 GeV lHOll . The nuclear modification factor Raa f° r inclusive J/xp from ALICE is almost 
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Figure 11: Invariant mass distributions of pairs from the same and mixed events in the 

J/ip region from ALICE [27]. The left and right panels represent the semi-central (50 - 80%) and 
the most central (0 - 10%) classes, respectively, in PbPb collisions at y / S]\nv = 2.76 TeV. 
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Figure 12: Nuclear modification factor Raa of inclusive J/ip as a function of N par t for pj > 
GeV/c and 2.5 < y < 4 in PbPb collisions at ^/s^ = 2.76 TeV from ALICE |13|. The statistical 
and systematic uncertainties are shown by bars and boxes, respectively. The ALICE data are 
compared to the PHENIX data (the open circles for midrapidity in |y| < 0.35 and the filled 
circles for the forward rapidity in 1.2 < |y| < 2.2) measured at low pj in AuAu collisions at 
= 200 GeV HU1. 



a factor of three larger than the PHENIX data at forward rapidity (1.2 < |y| < 2.2) for N part > 
180. In addition, differently from the PHENIX results, the ALICE data do not show centrality 
dependence. 



5 Comparison with Models 

There have been great efforts to understand the J/ip production in PbPb collisions at LHC. 
Figure 13 shows the comparison of the experimental data on Raa of J/tp with one of the trans- 
port model calculations as a function of centrality. This very transport model is based on the 
Boltzmann equation with the spectral properties of charmonia constrained by the thermal lat- 
tice QCD [28 J. The shadowing effect in the incoming nuclei has been introduced by artificially 



10 



5 Comparison with Models 
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Figure 13: Comparison of the experimental data on Raa of // tp with the transport model cal- 
culations by Zhao and Rapp as a function of N part in PbPb collisions at y/sj^ = 2.76 TeV l28l . 
The shadowing, Cronin, and formation- time effects are taken into account for the prompt J/tp 
production at midrapidity on the left panel. The recombined // xps are important at low pr for 
ALICE on the right panel, but negligible at high pj for CMS. 

lowering the cc cross section, da C c/ dy, in this (as well as other) transport model(s). As indicated 



in the left panel of Fig. 13 the shadowing effect somewhat reduces the prompt J/tp yields in 
the CMS acceptance (6.5 < pr < 30.0 GeV/c in |y| < 2.4), but the Cronin and formation-time 
effects are shown to enhance them. In the end, the model can reproduce the experimental data 
reasonably well except the most peripheral bin. On the other hand, recombination becomes the 
most important process, especially, at the low-p^ region. The right panel of Fig. 13 indicates 
that the recombination component can account for more than half of the inclusive J/tp yields 
in the ALICE acceptance (pr > GeV/ c in 2.5 < y < 4). It is noteworthy to point out that the 
ALICE data favor the existence of the shadowing effect in peripheral collisions, but disfavor it 
in central collisions. 



Figure 14 , taken from Ref . lH3ll , shows the comparison of the ALICE data with another transport 
model calculation that again employs the Boltzmann-type transport equation with the suppres- 
sion and regeneration parameters determined by the hydrodynamic equation |29|. (Note that 
the upper hatched area indicates the same range represented by the two solid lines in the right 
panel of Fig.[l3j) The lower and upper limits of each hatched area in Fig.[l4]represent with and 
without, respectively, shadowing. In both transport models, the regenerated J/tp component 
amounts about 50% of the measured inclusive // tp yield. 

In addition, the dashed lines in Fig.[l4]represent the calculations by the statistical hadronization 
model that assumes deconfinement and thermal equilibrium of the bulk of the cc pairs M30ll . In 
this model, charmonium can be formed at the phase boundary by statistical hadronization of c 
and c. In Fig.[l4j the calculations for two dc C c/ dy values (0.15 and 0.25 mb) are given. Similar 
with the transport models, it is found that the statistical hadronization model is also sensitive 
to the cc cross section in medium. 

Presently, a more detailed study for the cold nuclear matter effects on the J/tp production in 
heavy-ion collisions is underway using different parametrizations of the nuclear parton dis- 



tribution functions (nPDF). (Some examples of nPDF can be found in 1311 .) Figure 15 shows 
the comparison of the CMS and ALICE data on Raa of J/tp with the QCD calculations in the 
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Figure 14: Comparison of the ALICE data on Raa of inclusive / /\p with various model cal- 
culations as a function of N par t in PbPb collisions at s/s^ = 2.76 TeV HH. The hatched areas 
represent the two transport model calculations |28H29|. (The upper hatched area represents the 
same range as indicated by the two solid lines by Zhao and Rapp in the right panel of Fig. 13 ) 
The dashed lines represent the calculations by the statistical hadronization model 
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Figure 15: Comparison of the CMS (left) and ALICE (right) data on Raa of J/ip with the QCD 
calculations for the traditional 2—^2 kinematics at a parton level in PbPb collisions at v / snn = 
2.76 TeV. The blue (red) band was obtained by using the EKS (nDSg) parameterization for the 
modified nuclear parton distribution functions 11321 . 



traditional 2—^2 processes at a parton level (g + g —> J/xp + g) in PbPb collisions at a/s^n = 
2.76 TeV |32|. 

In addition, Fig.[l6]presents the comparison of the LHC data with the calculations obtained by 
the next-to-leading order (NLO) of the color evaporation model (CEM) before kj smearing l32ll . 
The overall cold nuclear matter effect caused by the modification of nPDF underestimates the 
suppression of the // ip yields relative to the pp data. (Only the nDSg parametrization in NLO 
of CEM barely touches the ALICE data at forward region.) The pPb run at LHC, being planned 



12 



6 Y from CMS 



5 

1.4 



1.2 
1 

0.8 
0.6 
0,4h 
0.2 



I 1 1 1 1 1 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I 1 1 1 1 I " 1 1 1 I 

CMS PbPb ys^ = 2.76 TeV 
■ Prompt J/y 




3 

a: 1.4 



CNM: CEM NLO 
I IEKS 

Lsnsberg. and Matayie [prelimirurp') ' 



-|y|<2.4 * + 

6.5 < p T < 30 GeV/c 

H i. I I I I I I I J ] I I [ L I I I I I I J I [ L L I I I J J I I I I I i L I I I I 

T) 50 100 150 200 250 300 350 400 



N. 



'cart 



- CNM; CEM NLO 
EKS 



i i i i I i i i i I i i i 



M | I I I I | I I I I | I I I I | i I I I 

PbPb \fs^ = 2.76 TeV 
* ALICE: J/y (2.5<y<4.0) J 

arXiv:1202.1383 




50 100 150 200 250 300 350 400 



N, 



p art 



Figure 16: Same as Fig. 15 but with the calculations obtained via the next-to-leading order of 
the color evaporation model l32ll . 

at the beginning of 2013, will be crucial to scrutinize the cold nuclear matter effects from hot 
final state effects in the heavy-ion data at LHC. 



6 Y from CMS 

Thanks to higher beam energy and luminosity, pp as well as heavy-ion collisions at LHC pro- 
vide a plenty of b quarks and, thus, bottomonia for the study of QGP. Although bottomonia are 
not the subjects of this workshop, we summarize the results on the Y family, as the physics is 
quite relevant to the charmonium production in heavy-ion collisions. 

Figure [T7| displays the nuclear modification factor Raa of Y(1S) as a function of N part in PbPb 
collisions at y/s^ = 2.76 TeV from CMS |15|. For the most central class (0 - 10%), R AA of Y(1S) 
is 0.45 ± 0.14(stat.) ± 0.08(syst.) ± 0.03(global) for < p T < 20.0 GeV/c, which indicates 
a suppression by about a factor two. However, large statistical and systematic uncertainties 
for less central event bins prevent drawing any definite conclusions on the centrality depen- 



dence. In the left panel of Fig. 17, the CMS data are compared to the preliminary STAR data for 
Y(1S+2S+3S) in AuAu collisions at y/s^j = 200 GeV ||33|. It is observed that a suppression of 
inclusive Ys at RHIC is similar to that at LHC despite the large gap in beam energies. 



In the right panel of Fig. 17 the CMS data are compared to the calculations by the anisotropic 
hydrodynamic model (AHYDRO) with different options for the ratio of the plasma shear vis- 
cosity to the entropy density (4717/ / S) |34|, where the term 'anisotropy' implies local momentum- 
space anisotropy due to finite plasma shear viscosity. The theoretical calculations for 4m]/S = 
1^3 can reproduce the measured centrality dependence of Raa f° r Y(1S). Note that, in the 
right panel of Fig. [17} the author of l34ll has assumed that about 51% of Y(1S) are produced 
directly and that all excited Y states have the same Raa as the Xbi- 



Figure 18 presents Raa of Y(1S) as functions of pr and y in PbPb collisions at ^/s^n = 2.76 TeV 
from CMS d. For the lowest p T bin from to 6.5 GeV/c, Raa of Y(1S) is 0.44 ± O.lO(stat) 
± 0.06(syst.) ± 0.04(global) in |y| < 2.4 for the minimum bias condition on centrality, which 
again indicates a suppression by about a factor of two, but the suppression disappears for 
Pt > 6.5 GeV/ c. For the rapidity dependence, a slightly larger suppression is observed at 
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Figure 17: Nuclear modification factor Raa of Y(1S) as a function of N part in PbPb collisions at 
yfs^j = 2.76 TeV from CMS [15J. The measured Ys are for < p T < 20 GeV/c in \y\ < 2.4. The 
grey box at Raa - 1 an d N part = represents a global scale uncertainty of 6% originating from 
the integrated luminosity of the pp data sample. The statistical and systematic uncertainties 
are shown by bars and boxes, respectively. (Left) The CMS data are in comparison with the 
Y(1S+2S+3S) data measured by STAR (|y| < 0.5) in AuAu collisions at ^/s^ = 200 GeV 1551 . 
The light-grey and open boxes at N part = 400 represent the systematic and statistical, respec- 
tively, uncertainties in pp collisions for STAR. (Right) The CMS data are compared with the 
theoretical calculations by the anisotropic hydrodynamic model with different options for the 
ratio of plasma shear viscosity to entropy density 4m] /S [34J. See the text for the details. 



midrapidity for minimum bias events, but the large statistical uncertainties do not allow us to 
make conclusions on the rapidity dependence at the moment. 



In Fig. 18 the CMS data are again compared with the calculations by AHYDRO with three 
different values of 4m] /S Il34ll . The calculations for 4mj/S = 1^3 can reproduce the measured 
rapidity dependence of Raa f° r Y(1S) reasonably well, but seem to be much flatter than the 
experimental data in the pj dependence. However, much more precise data are required to 
draw firm conclusions. 

Finally, the excellent resolution of the CMS detector for muons allows us to separate Y(1S) 
from the higher Y(2S+3S) states in heavy-ion collisions. Figure [19] shows the invariant mass 
distributions of ft + fi~ pairs for the pj of each muon, pj, larger than 4 GeV/c in pp as well 
as PbPb collisions at a/sjviv = 2.76 TeV lH6ll . Similar with the J/tp analysis, each Y state is 
parameterized by a 'crystal ball 7 function. Since the three Y states partially overlap in the 
measured mass range, a simultaneous fit with three 'crystal ball' functions added to the second- 



order polynomial function for the background has been performed. The solid lines in Fig. 19 
represent the final results of fitting. 

The comparison of the invariant mass distribution in pp with that in PbPb collisions 

clearly demonstartes that the higher Y(2S+3S) states are suppressed relative to the Y(1S) state 
in PbPb collisions. The fit functions to the pL + }i~ invariant mass distributions for pp and PbPb 



14 



7 Summary 



2.5 



a: 



1.5 



CMS PbPb ^ = 2. 
♦ V(1S) 



j | I I I | I I I | > I I 

76 TeV 



M. Strickland 
— 4nr|/S= 1 



0.5 



°0 



■ i I i i i I i 



I i i i I i i i I i i i I i 



Cent. 0-100% 

M<2.4 

I 1 I I I I I I I 



C£ 1.4 



1.2 
1 

0.8 
0.6 
0.4 
0.2 



' I M ' I ' ' ' I ' 



1 I " 1 I 1 



CMS PbPb \s m = 2,76 TeV 
♦ T(1S) 




8 10 12 14 16 18 20 

P T (GeV/c) 



- Cent 0-100% 
< p r < 20 GeV/c 



1 1 1 1 1 1 i 



uL: 



_ 



L_l 



M. Strickland 
— 4icn/S = 1 
— - 4Kt\fS ■ 2 

4jrf|/S - 3 

i I i i f 1 1 i i I i i 1 1 1 < 



2 0.4 0.6 8 1 1.2 1 4 1 6 1.8 2 2 2 2 4 



|y| 



Figure 18: Nuclear modification factor Raa of Y(1S), integrated over centrality, as functions of 
pr (left) and y (right) in PbPb collisions at y/s^jq = 2.76 TeV from CMS lH5t Raa is integrated 
over |y| < 2.4 for the pj dependence and over < pj < 20 GeV/c for the y dependence. The 
grey boxes at Raa = 1 represent global scale uncertainties of 6% originating from the integrated 
luminosity of the pp data sample. The statistical and systematic uncertainties are shown by 
bars and boxes, respectively. The CMS data are compared to the calculations by the anisotropic 
hydrodynamic model with different options for the ratio of plasma shear viscosity to entropy 
density Il34l . See the text for the details. 



give a double ratio 

Y(2S + 3S)/Y(lS)| PbPb n 

Y(2S + 3S)/Y(1S)| PP " 031 -o.i5(stat) ± 0.03(syst.), (4) 

where the systematic uncertainty is from varying the line shape in the simultaneous fit proce- 
dure, thus taking into account partial cancellation of systematic effects. 

In order to check the significance of Eq. |3J an ensemble of one million pseudo experiments 
has been generated with the signal line shape obtained from the pp data (the top panel of 
Fig. [19]), the background line shapes from both date sets, and a double ratio of Eq. [4] equal to 
unity within the statistical and systematic uncertainties. In this test, the probability to find the 
measured double ratio of 0.31 or less is estimated to be 0.9%, which corresponds to the 2.4 a 
effect in a one-tail integral of a Gaussian distribution. 

Currently, a great number of theoretical models are attempting to understand the Y production 
data from CMS. Some of them are successful reproducing the centrality dependence of Raa f° r 
Y(1S) and the double ratio given in Eq, gl I55H57I . 



7 Summary 

This article summarizes the experimental data on the J/tp production via channel in 

PbPb collisions at y/sjsfN = 2.76 TeV, measured by the ALICE, ATLAS, and CMS collaborations 
at LHC. Soon after the first heavy-ion run at LHC in 2010, ATLAS published the central-to- 
peripheral ratio of the J/tp yield, R cp , and demonstrated that more suppression occurs for more 
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(bottom) collisions at s/s^n = 2.76 TeV from CMS (16). The solid lines are the fit functions 
described in the text. 



central collisions. Later CMS published much more detailed results on the J/xp production us- 
ing the reference pp data. In the CMS analysis, the prompt and non-prompt // xps are separated 
by fitting the invariant mass and the pseudo-proper decay length distributions simulta- 

neously. The prompt and non-prompt J/xp yields were analyzed in PbPb collisions, and com- 
pared to those in pp collisions at the same beam energy in order to learn whether there exists 
any medium effect. 

The prompt // ip yield is suppressed by a factor of about 5 for the most central class (0 - 10%), 
which diminishes to about 1.6 for less central events (50 - 100%). The Raa of the prompt // ip 
at LHC is more than a factor two smaller than the J/xp data at RHIC. 

The Raa distribution of prompt J/xp displays no pj dependence with the current widths of 
the pr bins. In addition, slightly more suppression is observed at midrapidity than at forward 
rapidity, which is opposite to the behavior observed at low pj at RHIC. On the other hand, the 
non-prompt // xps decayed from the B mesons are also suppressed by a factor of about 2.7 with 
respect to the pp data. This result indicates the energy loss of fc-quarks in medium for the first 
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time. 

ALICE measured // ips at low pj in the forward region that is complimentary to the ATLAS and 
CMS acceptances. The low-pj Raa, integrated over centrality, shows a suppression relative to 
pp in the forward region. The comparison with the PHENIX data at forward rapidity indicates 
that R AA ofj/xpby ALICE is almost a factor of three larger than that by PHENIX for N par t > 180. 
Furthermore, the ALICE data do not show the centrality dependence contrary to the PHENIX 
data. 

The Raa distributions of prompt J/xp are compared with the calculations by the microscopic 
transport models and the statistical hadronization model. From these comparisons we learned 
that (1) regeneration is more important at low pj, (2) the theoretical calculations are sensitive 
to the cc cross section in medium, and (3) the calculations also depend on the nuclear parton 
distribution function in the incoming nuclei. Lastly, it has been demonstrated that the modifi- 
cation of nPDF alone cannot satisfactorily describe the LHC data on the J/xp yield. Presently, 
more detailed investigation on the cold nuclear matter effects in heavy-ion collisions is un- 
derway by testing different parametrizations of nPDF. The upcoming pPb run at LHC will be 
crucial to scrutinize the cold nuclear matter effects in the heavy-ion data. 

The Y(1S) state in PbPb collisions is suppressed by a factor of about 2 in the most central class (0 
- 10 %) when integrated over all pj and y. Suppression factor is large at low pj, but negligible 
at high pj. The ratio of Y(2S+3S) to Y(1S) in PbPb collisions is suppressed by a factor of about 3 
with respect to that in pp for pj > 4 GeV/ c. The probability to obtain the measured, or smaller, 
value of the double ratio, if the true suppression factor is unity, has been estimated to be less 
than 1%. Several theoretical efforts to understand the Y production data by CMS are ongoing. 
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